Scaling the spin densities separately in density functional theory 
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Coordinate scaling of each spin density separately is considered in spin density functional theory. 
A virial theorem relates the spin-scaled correlation energy to the spin-scaled correlation potentials. 
An adiabatic connection formula expresses energies at different spin interaction strengths in terms 
of spin scaling. Several popular approximate functionals are evaluated on the spin-scaled densities 
of atoms and of the uniform electron gas. The differences between this and uniform scaling are 
discussed. 
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I. INTRODUCTION 

Density functional theory combines accuracy and celer- 
ity in a computational scheme which has long been used 
by solid state physicists and has become popular in quan- 
tum chemistry p[. The Hohenberg-Kohn theorem j^] 
demonstrates that the electron density uniquely char- 
acterizes a ground state electronic structure problem so 
that the total energy is a functional of the density. This 
idea suggests the computationally expedient Kohn-Sham 
scheme || where only one part of the total energy func- 
tional must be approximated. This is the exchange- 
correlation energy, E xc , and the calculational accuracy 
of density functional theory is limited by the accuracy of 
approximations to -Exc! therefore, improvements in this 
functional are of great import. 

Exact constraints limit the possible forms of approxi- 
mations to _Bxc and provide guidance for the construc- 
tion of approximations Q. For total electron density 
functionals, Levy and Perdew || discovered a set of im- 
portant scaling and integral requirements that the exact 
functional must satisfy. They introduced the concept of 
uniform coordinate density scaling. This scaling takes a 
density n(r) into 



(r) = -f 3 n(jr), < 7 < 00, 



(1) 



and is a natural way to explore the behavior of density 
functionals. Many properties of the exact functional have 
been found by studying its behavior under this scaling 
H fj], ||. For example, the exchange energy changes as 



BxK] = lE x [n]. 



(2) 



All commonly-used approximations such as the local den- 
sity approximation (LDA), the generalized gradient ap- 
proximation (GGA), and hybrids of GGA with exact ex- 
change M satisfy this relation, some by construction. But 



for correlation, only inequalities can be derived, 

E c [n 1 ]> 1 E c [n] (7 > 1), (3) 

which some approximations satisfy |I| |l^| while others do 
not pl| . Levy and Perdew further showed that transla- 
tional invariance implies virial theorems relating poten- 
tials to energies. An example is 



E x [n] = - J d 3 r n(r) r • Vv x [n](r) 



(4) 



with a more complicated corresponding relation for cor- 
relation. The virials have been used to construct energy 
densities directly from potentials Exact statements 
about exact functionals are nontrivial and extremely use- 
ful in the construction and analysis of approximate func- 
tionals. Functionals that violate these exact conditions 
are unlikely to give reliable and physical results when 
applied to wide ranges of materials. 

On the other hand, modern density functional calcula- 
tions do not employ density functionals but rather use 
spin-density functionals. The basic idea is to replace 
-E xc [7i] with i? xc [nj, rtj so that the universal functional 
depends explicitly on both the up and down spin electron 
densities. Formal justification for this scheme was first 
given by von Barth and Hedin [|l3] and later by Rajagopal 
and Calloway |Q although some fundamental questions 
remain |l5| . There are several compelling reasons for 
using spin density functionals instead of total electron 
density functionals. Spin density functionals can more 
accurately describe systems with odd numbers of elec- 
trons [|l6|. They also allow the treatment of electrons 
in collinear magnetic fields and yield magnetic response 
properties p7| . Accurate calculation of these properties 
would be far more difficult in total-density functional the- 
ory, because a local spin density functional is a non-local 
total density functional. One cannot exaggerate how use- 
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ful spin density functional theory has been in accurately 
and efficiently calculating physical properties. 

To help develop improved spin-density functional, it 
would be of great interest to develop a formalism that 
probes the spin dependence of functionals and yields ex- 
act conditions about their spin-dependence. Here, we 
investigate whether scaling techniques developed for to- 
tal density functionals can be generalized to spin den- 
sity functional theory. Towards this end, we generalize 
uniform coordinate scaling of the density, Eq. ([j]) , to 
separate scaling of spin-densities: 

n-|a(r) = a 3 n|(ax), < a < oo 

n l/3 (r) = /3 3 n x (/3r), < p < oo. (5) 

According to this scheme, a spin-unpolarized system be- 
comes spin-polarized for a =/= (3. There are many other 
ways we could have chosen to alter the spin densities. 
For example, we could require that the total density re- 
main constant while the polarization changes. However, 
such a transformation would require introduction of en- 
sembles because fractions of electrons would be changing 
spin. The present scheme is simply the logical extension 
of coordinate scaling to separate spin densities. 

A principal result of this work is that, under such a 
transformation, a spin dependent virial theorem holds 
true. This theorem can be used to carefully check the 
convergence of spin-density functional calculations. It 
also provides a method for calculating the exact depen- 
dence of the correlation energy on spin-density for model 
systems for which accurate Kohn-Sham potentials have 
been found. Finally, it can be integrated over the scale 
factor a to give a new formal expression for the func- 
tional, E xc . 

Considerable progress has been made in density func- 
tional theory by writing E xc as an integral over a cou- 
pling constant A in what is called the adiabatic connec- 
tion relationship |l8|, |l9). For example, the success of 
hybrid functionals such as B3LYP |)[ ^(J can be under- 
stood in terms of this adiabatic connection pi], ^] . The 
adiabatic connection is simply related to uniform coordi- 
nate scaling |2^, ^4|. By analogy, we relate spin scaling 
to a spin-coupling constant integration, and we define a 
suitable generalization for this definition with a coupling 
constant for each spin density. 

We illustrate our formal results with several cases. For 
the uniform electron gas, we can perform this scaling 
essentially exactly. We show how this transformation re- 
lates energies to changes in spin-polarization. In this 
case, considerable care must be taken to deal with the 
extended nature of the system. We also show the results 
of spin scaling densities of small atoms using presently 
popular approximations. We close with a discussion of 
the fundamental difficulty underlying this spin scaling 
approach. 

Throughout, we use atomic units (e 2 = Ti = m e = 1), 
so that all energies are in Hartrees and all lengths in Bohr 
radii. We demonstrate all scaling relationships by scaling 
the up spin densities. Results for scaling the down spin 



are obtained in a similar fashion. Just change the spin 
label! 



II. SEPARATE SPIN SCALING THEORY 

The first interesting property of the spin scaling trans- 
formation, Eq. (|^), is that it conserves the total number 
of electrons globally even though the scaled spin density 
might tend towards zero at any point. As a diminishes, 
the two spin densities occupy the same coordinate space, 
but on two very distinct length scales. Even when a — > 0, 
the up electrons do not vanish, but are merely spread over 
a very large volume. The scaled density presumably then 
has vanishingly small contribution to the correlation en- 
ergy. For finite systems, we can consider this limit as the 
effective removal of one spin density to infinitely far way. 
We will discuss what this means for extended systems 
later when we treat the uniform gas. 

Another interesting property of the spin scaling trans- 
formation is that a scaling of one spin density can always 
be written as a total density scaling plus an inverse spin 
scaling of the other spin; that is 

E xc [n 1a ,n L ] = E xc [{ n h n ll/a } a ] (6) 

where the parenthesis notation on the right indicates scal- 
ing the total density. Thus, without loss of generality, we 
need only scale one spin density. 

To understand what happens when a single spin den- 
sity is scaled, we first study exchange. Because the spin 
up and down Kohn-Sham orbitals are independent, the 
exchange energy functional can be split into two parts, 
one for each spin |^5| . The scaling relationships for total 
density functional theory generalize for each term inde- 
pendently. For an up spin scaling, we find 

= |l? x [2n T ] + \Ex\2h\. (7) 

When a — > 0, we are left with only the down contribution 
to exchange. Separate spin scaling allows us to extract 
the contribution from each spin density separately, e.g., 
dE x [n^ a , nj/da at a = 1 is the contribution to the ex- 
change energy from the up density. A plot of E x [n-\ a) nj 
versus a between and 1 yields a straight line and is 
twice as negative at 1 as at 0. 

Separate spin-scaling of the correlation energy is more 
complicated. Unlike E x [n^, nj, E c [n^ , n J cannot triv- 
ially be split into up and down parts. The Levy method 
of scaling the exact ground-state wave-function does not 
yield an inequality such as Eq. (j^) because the spin- 
scaled wave-function is not a ground-state of another 
Coulomb-interacting Hamiltonian. Nor does it yield an 
equality as in the spin-scaled exchange case, Eq. (Q), 
because the many-body wave-function is not simply the 
product of two single spin wave- functions. In both cases, 
the two spins are coupled by a term l/|r — ar'\. 
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To obtain an exact spin scaling relationship for E xc , 
we take a different route. Consider a change in the energy 
due a small change in the up-spin density: 

SE XC = E xc [n^ + Sni,ni\ - E xc [n h n{\. (8) 

Use t>xcT( r ) = 5E XC /Sni(r) to rewrite 5E XC as 



5E XC = / d 3 r <5n T (r) v xcT [n T ,r«x](r). 



(9) 



to first order in Sn-\. Now, consider this change as com- 
ing from the following scaling of the density, n^ a (r) — 
a 3 n|(ar), where a is arbitrarily close to one. The change 
in the density is related to the derivative of this scaled 
density: 



dn ]a (v) 
da 



a =i = 3n-f (r) + r • Vn T (r) . 



(10) 



Use Eqs. (|9|) and (|10|), and integrate by parts to find 

dE xc [n Ta ,n|] 



da 



■|o=i = - J d 3 r ?i T (r) r • Vv xcT [n T , raj(r). 



(11) 

Eq. ( |ll| ) is an exact result showing how dE xc /da\ a= i 
can be extracted from the spin densities and potentials. 
For an initially unpolarized system, rz.f = = n/2, and 
^xcT = v xci — Vxc ■ Thus the right-hand-side of Eq. (|ll| ) 
becomes half the usual virial of the exchange-correlation 
potential. This virial is equal to dE xc [n a ]/da\ a= i = 
£ xc + T c . T c is the kinetic contribution to the corre- 
lation energy^). Thus, for spin-unpolarized systems, 



dE c [n Ta ,nJ 
da 



1 



(E c + T c ) 



(12) 



For initially polarized systems, there is no simple relation 
between the two types of scaling. 

To generalize Eq. ( |ll| ) to finite scalings, simply replace 
ri| on both sides by n^ a , yielding: 

dE xc [n^ a ,n{\ 1 f 3 

= / d°r n ]a {r) r • Vv xcT [n Tce , n x J (r). 

da a J 

(13) 

We can then write the original spin-density functional as 
a scaling integral over this derivative: 

t? r l i- j? r l , /"^ ^xcI^Ta^l] 

A xc [n T ,nJ = hm ii c [?i T „,nJ + / da- 



(14) 

This is a new expression for the exchange-correlation en- 
ergy as an integral over separately spin-scaled densities 
where the spin-scaled density is scaled to the low-density 
limit. With some physically reasonable assumptions, we 
expect 



lim £ c [n TQ ,nJ = £ c [0,nj. 



(15) 



For example, if the anti-parallel correlation hole vanishes 
as rapidly with scale factor as the parallel-spin correla- 
tion hole of the scaled density, this result would be true. 



Numerical results indicate that this is the case for the ap- 
proximate functionals used in this paper. Nevertheless, 
we have not proven Eq. ( |l5| ) here. 

A symmetric formula can be written down by scaling 
the up and down spins separately and averaging: 

E xc K,nJ = - lim (£; xc [n TQ ,nJ + E xc [n h n la }) 

Z a— »0 

+ — / da. I d?r n^ a (r)r • Vi? X cT [ n 1a ? 

+ \J o d(5 J d 3 r n l0 (r)r ■ Vw xc| [ti t , n i/3 ](r)(16) 

This result is the spin density functional generalization 
of spin-decomposition, coordinate scaling, and the virial 
theorem. Each of these ideas yields separate results for 
pure exchange or uniform coordinate scaling, but all are 
combined here. Notice that the potentials depend on 
both spins, one scaled and the other unsealed. This re- 
flects the difficulty in separating up and down spin cor- 
relations. 

The proof of Eq. (0) is true for exchange-correlation, 
but in taking the weakly-correlated limit, the result also 
holds true for exchange. In the exchange case, Eq. ( |l6| ) 
reduces to Eq. (|?]) with equal contributions from the limit 
terms and the virial contributions. To obtain this result, 
recall how E x scales, Eq. (||). Since the energy contri- 
bution from each spin is separate and since the sca ling 
law is linear, the limits in the first two terms of Eq. (|lq ) 
are doable without any extra physical assumptions. The 
virial terms are a bit more difficult to handle as the ex- 
change potentials change under scaling. In the end, the 
first two terms contribute half the exchange energy while 
the virial terms contribute the other half. 



III. UNIFORM GAS 

To illustrate the effect of spin scaling, we examine the 
uniform electron g system for which we have essen- 
tially exact results. Great care must be taken to define 
quantities during separate spin-scaling of extended sys- 
tems. Begin with a spin-unpolarized uniform electron 
gas of density n and Wigner-Seitz radius r s = (3/47rn)3 . 
When one spin density is scaled, the system becomes 
spin-polarized, and relative spin-polarization is measured 
by 



C 



(17) 



We assume that for a spin-polarized uniform system, the 
exchange-correlation energy per electron, e x " (r s ,0, is 
known exactly. We use the correlation energy parame- 
terization of Perdew and Wang |Q to make our figures. 

To perform separate spin scaling of this system, we 
focus on a region deep in the interior of any finite but 
large sample. A simple example is a jellium sphere of 
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FIG. 1: Spin scaling of a uniform gas: exchange energy per 
particle, ex(?if Q , nj.), at r a — 2 (dotted line) and 6 (solid line). 
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FIG. 2: Spin scaling of a uniform gas: correlation energy per 
particle, e c (n^ a ,ni), at r s = 2 (dotted line) and 6 (solid line). 



radius R » r s . The correlation energy density deep in 
the interior will tend to that of the truly translationally 
invariant uniform gas as R — ► oo. At a = 1, we have an 
unpolarized system with n-p = = n/2. The up-spin 
scaling, n-f = a 3 n/2, changes both the total density and 
the spin-polarization. Deep in the interior 



r s (a) 



1 



1/3 



(18) 



where r s is the Seitz radius of the original unpolarized 
gas, and 



a 6 -1 



(19) 



ct 3 + 1' 

The energy density here is then 

exc(a) = elf (n Ta , nj = e£*(r.(a), ((a)), (20) 



and the energy per particle is 

e X c(ai) = e X c(aO/n(a) 



(21) 



where n(a) is the interior density. 

To illustrate the effects of this spin scaling, consider 
the simple exchange case. Deep in the interior, we have 
a uniform gas of density n-\ a and rij, and the energy den- 
sities of these two are given by Eq. (^) , since the integrals 
provide simple volume factors. The Slater factor of n 4 / 3 
in the exchange density of the uniform gas produces a 
factor of (1 + a 4 ). When transforming to the energy per 
electron, there is another factor of (1 + a 3 ) due to the 
density out front. Thus the exchange energy per electron 
is 



ex(a) 



1 



1 



unpol. 



(n) 



(22) 



This variation is shown in Fig. |l|. This result may ap- 
pear to disagree with Eq. ([?]), but it is valid deep in the 
interior only. To recover the total exchange energy, one 
must include those electrons in a shell between R and 
R/a with the full polarized uniform density a 3 n/2. The 
exchange energy integral includes this contribution, and 
then agrees with Eq^Jffl). 

Near a = 1, Eq. (§|) yields (1 + a) eT P ° L A in agree- 
ment with a naive application of Eq. (Q) . This is because, 
in the construction of the energy from the energy per 
electron, the factor of the density accounts for changes 
in the number electrons to first order. So the derivative 
at a = 1 remains a good measure of the contribution to 
the total exchange energy from one spin density. On the 
other hand, as a — > 0, the exchange energy per electron 
in the interior returns to that of the original unpolarized 
case. This reflects the fact that exchange applies to each 
spin separately, so that the exchange per electron of the 
down-spin density is independent of the presence of the 
up-spin density. 

Figure || shows the uniform electron gas correlation 
energy per particle scaled from unpolarized (a = 1) to 
fully polarized limits (a = 0). Again, the curves become 
flat as a — > 0, because for small a, there is very little 
contribution from the up-spins. Now, however, there is a 
dramatic reduction from a — 1 to a — because of the 
difference in correlation between unpolarized and fully 
polarized gases. Note that the correlation changes tend 
to cancel the exchange variations. 



IV. FINITE SYSTEMS 

Next, we examine the behavior of finite systems under 
separate-spin scaling. We choose the He and Li atoms, 
to demonstrate the effects on the simplest non-trivial un- 
polarized and spin-polarized cases. For each system, we 
solve the Kohn-Sham equations using a specific density 
functional approximation. The resulting self-consistent 
densities are then spin scaled and the approximate en- 
ergies evaluated on the scaled densities using that same 
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FIG. 3: Spin scaling of the He atom density using vari- 
ous approximate functionals for Ec'- local spin density ap- 
proximation (solid line), generalized gradient approximation 
(PBE, dashed line), BLYP (bars), self-interaction corrected 
LSD (short dashes). 



functional. Since these are approximate functionals, nei- 
ther the densities nor the energies are exact. We are 
unaware of any system, aside from the uniform gas, for 
which exact spin-scaled plots are easily obtainable. For 
now, we must compare plots generated from approxi- 
mate functionals. Even the simple atomic calculations 
presented here were rather demanding since, especially 
for very small spin-scaling parameters, integrals contain- 
ing densities on two extremely distinct length scales are 
needed. 

The He atom (Fig. |^) is spin unpolarized at a = 1. 
Scaling either spin density gives the same results. The 
LSD curve gives far too much correlation and does not 
vanish as a — ► 0. In the fully scaled limit, we expect, 
as we have argued in section El that the correlation en- 
ergy should vanish. This is because the two electrons 
are now on very different length scales and so should not 
interact with each other. The residual value at a — > 
reflects the self-interaction error in LSD for the remain- 
ing (unsealed) one-electron density. The PBE curve is 
on the right scale but also has a residual self-interaction 
error as a — > 0. The BLYP functional 11, |6) is pop- 
ular in quantum chemistry and gets both limits correct. 
However, the functional's lack of self-interaction error is 
because the correlation energy vanishes for any fully po- 
larized system. This vanishing is incorrect for any atom 
other than H or He. Finally, the LSD-SIC curve 
is probably the most accurate in shape (if not quantita- 
tively) since this functional handles the self-interaction 
error appropriately. We further observe that the curves 
appear quite different from those of the uniform gas. The 
atomic curves are much flatter near a — > 1 and have ap- 
preciable slope near a — > 0. This is because these ener- 
gies are integrated over the entire system, including the 



Approx. 


E x 


E c 


E c [n r ,0\ 


dE c /da 


LSD 


-0.862 


-0.111 


-0.027 


-0.022 


PBE 


-1.005 


-0.041 


-0.005 


-0.002 


SIC 


-1.031 


-0.058 


0.000 


-0.011 


BLYP 


-1.018 


-0.044 


0.000 


-0.005 


exact 


-1.026 


-0.042 


0.000 


-0.003 



TABLE I: He atom energies, both exactly and within sev- 
eral approximations. All energies in Hartrees; all functionals 
evaluated on self-consistent densities. 
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FIG. 4: Up spin scaling of the Li atom density using various 
approximate functionals for Ec'- local spin density approxi- 
mation (solid line) , generalized gradient approximation (PBE, 
dashed line), BLYP (bars). 



contribution from the entire spin-scaled density, whereas 
the energy densities in the uniform gas case were only 
those in the interior. 

Quantitative results are listed in Table Q. The exact He 
values, including the derivative at a = 1, using Eq. jl^), 
were taken from Ref. @ Note that PBE yields the 
most accurate value for this derivative. The BLYP cor- 
relation energy is too flat as function of scale parameter. 
BLYP produces too small a value for T c leading to a lack 
of cancellation with E c and a subsequent overestimate 
of the derivative at a = 1. LSD-SIC has a similar prob- 
lem. The LSD value, while far too large, is about 8% of 
the LSD correlation energy, close to the same fraction for 
PBE, and not far from exact. However, the important 
point here is that results from separate spin-scaling are a 
new tool for examining the accuracy of the treatment of 
spin-dependence in approximate spin-density functionals. 

The Li atom (Figsjj and ||) is the smallest non-trivial 
odd-electron atom. We choose the up spin density to 
have occupation ls2s. As the up spin is scaled away, as 
in Fig. E, we find a curve very similar to that of He, 
Fig. pL The primary difference is the greater correlation 
energy for a = 1. 

On the other hand, scaling away the down-density 
gives a very different picture, Fig. 0. The most dramatic 
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FIG. 5: Down spin scaling of the Li atom density using var- 
ious approximate functionals for Ec'- local spin density ap- 
proximation (solid line), generalized gradient approximation 
(PBE, dashed line), BLYP (bars). 



Approx. 


E x 


Eg 


E c [n u 0] dEc/da 


E a [0, ni ] 


dEc/d/3 


LSD 


-1.514 


-0.150 


-0.047 -0.037 


-0.032 


-0.019 


PBE 


-1.751 


-0.051 


-0.012 -0.004 


-0.005 


-0.001 


BLYP 


-1.771 


-0.054 


-0.054 -0.020 


0.000 


-0.005 


exact 


-1.781 


-0.046 




0.000 





TABLE II: Li atom energies, both exactly and within sev- 
eral approximations. All energies in Hartrees, all functionals 
evaluated on self-consistent densities. 



changes in the correlation energy now occur at small a. 
Near a — ► 1, the system energy is quite insensitive to 
spin-scaling, especially in GGA. This is exactly oppo- 
site to what we have seen for the uniform gas. It is an 
open question whether this would be observed with the 
exact functional. For up spin scaling, we expect the cor- 
relation energy to vanish as a — * 0. But for down spin 
scaling one expects a finite correlation energy in the limit 
/3 — * 0. The two spin- up electrons remain and are still 
correlated. In this case, the BLYP functional errs no- 
ticeably since it predicts no correlation energy for the 
remaining two electrons. 

Quantitative results for Li are given in Table H The 
exact result for E x is the E x of a self consistent OEP 
calculation. Using the highly accurate energy prediction 
from |3^|, we deduce the exact E c = Et — Et.oep- The 
other exact results are not extractable from the literature 
here, but could be calculated from known exact poten- 
tials and densities |3lj]. Even in this simple case, an SIC 
calculation is difficult. For the up-spin density, one would 
need to find the Is and 2s orbitals for each value of a 
that yield the spin-scaled densities. 



V. SPIN ADIABATIC CONNECTION 

Here, we define an analog of the adiabatic connec- 
tion within the spin-scaling formalism. Traditionally, 
we think of A as a parameter in the Hamiltonian, but 
this way of thinking becomes prohibitively complicated 
in spin density functional theory. We would have to de- 
fine three coupling constants: Af, Aj,, and Afj, . Even if we 
did that, it would be non-trivial to relate changes in these 
coupling constants to changes in the electron density. In- 
stead, we define a relationship between spin-scaling and 
a spin dependent coupling parameter. For total density 
scaling, the relationship between scaling and evaluating 
a functional at a different coupling constant is P, p3| 



E$ [n]=\*E xa [nyx] 
The adiabatic connection formula is 



dX 



dEi 



dX 



d\ U XC (X). 



(23) 



(24) 



By virtue of the Hellmann-Feynman theorem, U XC (X) can 
be identified as the potential contribution to exchange- 
correlation at coupling constant A. The integrand U XC {X) 
can be plotted both exactly and within density functional 
approximations, and its behavior lends insight into defi- 
ciencies of functionals |32) . For separate spin scaling, we 
apply the same ideas but now to 



A£ xc [n T ,nJ = E xa [n h ni\ - E xc [0,n;], 



(25) 



the exchange-correlation energy difference between the 
physical system and the system with one spin density 
removed while keeping the remaining spin-density fixed. 
For polarized systems, this quantity depends on which 
spin density is removed. We define 



AE x l = A?A£ X 



>[ n i/x v n l\ 



and 



so that 



AC/ XC (A T ) = dAE x l/dX h 



AE XC = / dA T AC/ XC (A T ). 
Jo 



(26) 



(27) 



(28) 



This produces a spin-dependent decomposition of the 
exchange-correlation energy, related to separate spin- 
scaling rather than total density scaling, with the integral 
now including the high-density limit. As Af — > 0, ex- 
change dominates, and f7 xc (Af ) — + U x (X-\ ) which is just 
E x [2ni]/2 according to the simple results for exchange 
in Sec. I. Furthermore, in the absence of correlation, 
t/ xc (A|) is independent of Af. This is not true if one uses 
a naive generalization of Eq. ( p^ ) . 

This spin adiabatic connection formula should prove 
useful for the improvement of present-day functionals in 
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FIG. 6: Single spin adiabatic connection for He atom: 
local spin density approximation (solid line), generalized 
gradient approximation (PBE, dashed line), BLYP (bars), 
self-interaction corrected LSD (short dashes), Exact (fancy 
dashes) . 





LSD 


PBE 


SIC 


BLYP 


exact 


Af/xc(0) 
AKtc(l) 

AiJxc 


-0.43 
-0.58 
-0.54 


-0.50 
-0.57 
-0.54 


-0.53 
-0.62 
-0.57 


-0.51 
-0.59 
-0.55 


-0.51 
-0.60 
-0.56 



TABLE III: Spin adiabatic connection A[/ X c(Af) for He 
atom, both exactly and in several approximations. 



the same way that the adiabatic connection formula has 
been useful for improving total density functionals. For 
example, it might be possible to perform Gorling-Levy 
perturbation theory ]33[ ] in this parameter (Af) or to ex- 
tract a correlation contribution to kinetic energy |3^| . 
We show the spin adiabatic connection for the He atom 
in Fig. ^. In generating each adiabatic connection plot, 
we now take the scaled spin density to the high density 
limit. The area under each curve is precisely AE XC for 
a particular approximation. To get E xc [ni,ni\, we must 
add the contribution from the unsealed spin, Z? xc [0,nj. 
The spin adiabatic connection curve looks quite similar to 
the regular adiabatic connection curve: for the He atom, 
AU XC (X^) becomes more negative with A everywhere and 
is close to linear. This suggests that the spin-correlation 
effects are weak for this system, just as the correlation 
effects are. 

To better understand how popular approximations 
perform, we would like to compare with the exact curve. 
In principle, this requires a sophisticated wavefunction 
calculation designed to reproduce the spin-scaled den- 
sities at every point in the adiabatic connection curve. 
Here, we use a simple interpolation that should be 
highly accurate. Analytic formulae give exact limits for 
AU K c(X-\)- At the small Af limit, exchange dominates, 
and we are left with the exchange contribution from the 



scaled spin to the total energy: 

A[/ XC (A T =0) = i^xM (29) 

At the other end, 

AC/ XC (A T = 1) = 2E xc [n h n L ] - 2£ xc [0,n ; ] 

-dE X c[ni a ,ni\/da\ a =i (30) 

For a spin-unpolarized two electron system like the He 
atom, this becomes 



A£/ XC (A T = 1) = E x /2 + 2E C - {E c + T G )/2 
(2 electrons, unpol.) 



(31) 



For He at Af = 1, AU XC (X] ) = —0.60. To approximate 
the exact curve, we use a (1,1) Pade approximant. The 
values AU xc (0), AU XC (1), and Ai5 X c nx the three un- 
known parameters. This pade turns out to be nearly a 
straight line. 

Table III shows the exact limits and the limits given 
by several popular functionals. BLYP reproduces both 
limits most accurately and is mostly linear. This should 
come as no surprise as BLYP yields good energies and 
accounts for He's self-interaction error (if a bit serendip- 
itously). However, we do not expect such good results 
from BLYP when using it on Li. For Li, as we have see 
in section IV , BLYP predicts no correlation energy when 
only one electron is scaled away. BLYP will fail notice- 
ably and uncontrollably in this case. The LSD functional 
dramatically underestimates the single spin exchange en- 
ergy and, therefore, gets the small Af limit quite wrong. 
This reflects the usual error for LDA exchange. But no- 
tice how well LSD performs performs at Af = 1. The 
value here is only a 3% overestimate of the exact value, 
much better than the 9% overestimate for the exchange- 
correlation energy. Furthermore, the LSD derivative as 
Aj — > 1 is almost exact. PBE and LDA-SIC are qualita- 
tively similar, the greater error in LSD-SIC being due to 
the errors in LSD. Both show a flattening of the curve as 



A 



0. much more than BLYP. Our exact curve is too 



crudely constructed to indicate which behavior is more 
accurate. 

Ideally, we would compare approximations to the exact 
adiabatic plot for this and other systems such as the Li 
atom. The plots are not easy to generate. But even so, 
analysis of the exact limits is sufficient to garner a deeper 
understanding of how functionals treat and mistreat spin 
densities. 



VI. CONCLUSIONS 

Both scaling and the adiabatic decomposition formula 
have proven extremely useful in studying and construct- 
ing total density functionals. We have suggested the pos- 
sibility of scaling spin densities separately, derived a new 
virial theorem, given new exact results for the He atom, 
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and pointed out the difficulties of deducing exact theo- 
rems from this decomposition. While exact calculations 
are difficult to perform and and exact results appear dif- 
ficult to prove within this approach, any results would be 
very useful and likely to improve spin density functional 
theory's treatment of magnetic properties. 

We close with a significant challenge to developing sep- 
arate spin scaling. In the total density scaling of Eq. 
([!]), the density is both squeezed (or spread) and is also 
translated. The squeezing is independent of the choice of 
origin, but the translation is not. This origin-dependence 
should not affect the exchange-correlation energy because 
space is translationally invariant. However, when an in- 
dividual spin density is scaled, the remaining spin density 
remains fixed in space. This means the resulting density 
depends on the choice of origin for the separate spin- 
scaling. So while E c [n-\ a , nj is a spin-density functional 
of n-\ a and rt^, it is not a pure spin-density functional 
of the original spin-densities because of this origin de- 
pendence. Most likely, a method of transforming away 
this origin dependence, as found for virial energy densi- 
ties in Ref . pj| , will be needed to make this spin scaling 
technique more physical and useful. For atoms, we made 
the obvious choice of origin at the center of the nucleus. 
Origin dependence will become acute in applications to 
molecules and even worse for solids. On the other hand, 
the non-uniform coordinate scaling of Gorling and Levy 



pq| suffers from the same difficulties for non-spherical 
densities but has still produced useful limits for approx- 



imate density functionals 37 



However, it is important to stress that the spin virial 
relationship is unaffected by this challenge. For a arbi- 
trarily close to 1, the spin-scaled energies are independent 
of the choice of origin, and these difficulties are irrelevant. 
The spin virial relationship is an exact constraint and 
give us a useful measure of how the correlation energy is 
affected by small changes in the spin densities. It also 
leads to a natural decomposition of energy changes due 
to separate spin densities. It should be useful in deter- 
mining whether calculations are self-consistent for each 
spin density separately. This might be useful for example 
in systems where small differences between spin densities 
are important to calculate properly. 
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